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In this paper we investigate the possibility of measuring the post-Newtonian general 
relativistic gravitomagnetic Lense-Thirring effect in the Jovian system of its Galilean 
satellites lo, Europa, Ganymede and Callisto in view of recent developments in processing 
and modelling their optical observations spanning a large time interval (125 years). The 
present day best observations have an accuracy between several kilometers to few tens of 
kilometers, which is just the order of magnitude of the Lense-Thirring shifts of the orbits 
of the Galilean satellites over almost a century. From a comparison between analytical 
development and numerical integration it turns out that, unfortunately, most of the 
secular component of the gravitomagnetic signature is removed in the process of fitting 
the initial conditions. Indeed, an estimation of the magnitude of the Lense-Thirring effect 
in the ephemerides residuals is given; the resulting residuals have a maximum magnitude 
of 20 meters only (over 125 years). 
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1. Introduction 

1.1. The post-Newtonian gravitomagnetic Lense-Thirring effect 

One of the few predictions of the Einstein's General Theory of Relativity (GTR) 
which is still awaiting for a direct and reliable observational check is the so called 
Lense-Thirring (LT) effect^. It is a consequence of the off-diagonal components 
goi, i = 1, 2, 3 of the space-time metric generated by a weakly gravitating and slowly 
rotating massive body of mass M and proper angular momentum S in the linearized 
weak- field and slow- motion approximation of GTR'^'"^. The mass-energy currents 
induce, among other things, a Lorentz-like, velocity-dependent force which acts on 
the orbital motion of a test particle freely moving in the gravitational field of the 
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central body. Such force is 

F=-2m(^) xBg, (1) 

where v is the particle's velocity and 

^ G[3rir.S)-r^S]_ 

cr° 

G is the Newtonian gravitational constant and c is the speed of light in vacuum. Due 
to the formal resemblance of Bg with the dipolar magnetic field of the Maxwellian 
electromagnetism, the ensemble of gravitational phenomena induced by the mass 
currents is named gravitomagnctism^ . It turns out that the non-central gravitomag- 
netic force of Eq. (1) makes the longitude of the ascending node and the argument 
of pericentre w to undergo tiny secular precessions^'*"^ 

dVLuT 2GS diJuT 6GScosi 



dt c2a3(i_e2)3/2' dt c2a3(l-e2)3/2' 



(3) 



where a, e, i are the semimajor axis, the eccentricity and the inclination yto the 
equator of the orbit of the moving particle, respectively. 

The gravitomagnetic phenomena may have strong consequences in many astro- 
physical and astronomical scenarios involving, e.g., accreting disks around black 
holes^, gravitational lensing and time delay^'^. Unfortunately, in these contexts the 
knowledge of the various competing effects is rather poor and makes very difficult 
to reliably extract the genuine gravitomagnetic signal from the noisy background. 
E.g., attempts to measure the LT effect around black holes are often confounded by 
the complexities of the dynamics of the hot gas in their accretion disks. Conversely, 
if one believes in GTR, the LT effect could be used, in principle, for measuring 
the proper angular momentum of the central mass which acts as sources of the 
gravitomagnetic field. 

It must be noted that, according to K. Nordtvcdt, gravitomagnetism would have 
already been indirectly tested from the radial motion of the Earth's geodetic satellite 
LAGEOS^*^ and from the high-precision reconstruction of the lunar orbit with the 
Lunar Laser Ranging (LLR) tcchnique^^. 

An attempt to directly measure the LT effect in the gravitational field of the 
Earth has been performed by CiufoUni and PavUs with the LAGEOS and LAGEOS 
II satellites'^ by using an observable proposed by one of us in Ref. 13. The claimed 
total accuracy is 5-10% at 1-3 sigma, but such estimates have been criticized by 
one of us in Ref. 14 who proposes a 19-24% total error at 1-sigma. The main 
limiting factors are the systematic errors of gravitational origin due to the static and 
secularly varying parts of the even zonal coefficients of the multipolar expansion 
of the terrestrial gravitational potential. 

The first, very preliminary, evidence of the gravitomagnetic field of the Sun 
through the LT precessions of the longitudes of the perihelia w of the inner planets 
of the Solar System has recently been reported by one if us in Ref. 15. 
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In April 2004 the GP-B spacecrafts^ has been launched. Its aim is the measure- 
ment of another gravitomagnetic effect, i.e. the precession of the spins^^ of four 
superconducting gyroscopes carried onboard with an expected accuracy of 1% or 
better. 



1.2. Aim of the paper 

In this paper we want to investigate the possibility of measuring the LT effect on 
the orbits of the Galilean satellites lo (1), Europa (2), Ganymede (3) and Callisto 

(4) evolving in the gravitational field of Jupiter in view of recent improvements in 
their ephemeridcs^^'^^. The orbital parameters of the Galilean satellites of Jupiter 
are presented in Table 1. It is noteworthy that the idea of using the Jovian system 

Table 1. Averaged Kcplcrian orbital ele- 
ments of the Galilean satellites lo, Europa, 
Ganymede and Callisto deduced from LI 
ephemerides. 



Satellite 


a (km) 


e 


i (dog) 


lo 


422,030 


0.0042 


0.036 


Europa 


671,261 


0.0094 


0.469 


Ganymede 


1,070,621 


0.0015 


0.175 


Callisto 


1,883,134 


0.0075 


0.187 



was put forth for the first time by Lense and Thirring themselves in their original 
paper s. 



2. The R-T-N scheme for the Lense-Thirring effect 

Since the Galilean satellites of Jupiter move along nearly circular and equatorial 
orbits and to make easier the comparison with the current ephemerides we will 
consider the radial, transverse and normal components of their orbits. Their per- 
turbations can be expressed in terms of the integrated shifts of the Keplerian orbital 
elements Aa, Ae, Ai, Afl, Aw, AAA, where A4 is the mean anomaly, as^*^ 



AR = \ {Aaf + 



[(eAa + aAe)2 + [aeAMf] 



AT = a\ 1 + 



AN : 



AM + Aw + cosiAfl + ^/ {Aef + {eAMf 




{Aty 



+ {siniAny^ 



(4) 
(5) 
(6) 



^Prom Eq. (3) it turns out that the gravitomagnetic force only affects, in general, 
both the transverse and the normal components: Ai?LT = 0. For polar orbital 
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configurations, i.e. i = 90 deg, ATlt = 0, while for equatorial orbits, i.e. i = deg, 
AJVlt = 0. For nearly circular and equatorial orbits it turns out that AT ~ aAA, 
where X = M + fl + co is the mean longitude. 

3. The Lense-Thirring effect in the Jovian system 

The four Galilean satellites'' of Jupiter arc currently the best candidates (among 
satellites evolving around giant planets) in regard to their orbital motion's accuracy. 
Indeed, ephemerides of the Galilean satellites reach an accuracy of only few tens 
of kilometers for the modern period^^, and benefit of observations dispatched over 
more than a century (from 1891 to 2003 in the case of LI ephemerides*^). These 
observations have various origin. While most of them are photographic ones, the 
observation of mutual events (the eclipse or occultation of one satellite by another) 
delivers the most accurate Earth-based observations of the Galilean system. Indeed, 
such observations are photometric instead of astrometric, and so are less sensitive 
to atmosphere turbulence. Hence, the accuracy of these latter can reach few tens of 
kilometers. 

More, JPL ephemerides accuracy amount to 5 km for the present pe- 
riod by the use of unpublished spacecraft observations (see Ref. 21 and 
http: / / ssd.jpl.nasa.gov/sat_eph.html) . 

3.1. Theoretical predictions: analytical and numerical estimates 

In regard to the LT effect, for the Jupiter's angular momentum we assume 
S = 4.33 X 10^^ kg m^ s~^ since the ratio'' a of the moment of iner- 
tia I to MR"^, where R is the mean equatorial radius, amounts to^^ 0.264 
(see also http://nssdc.gsfc.nasa.gov/planetary/factsheet/jupiterfact.html). Such es- 
timates are based on theoretical models: at present, there are no direct, independent 
measurements of S for Jupiter. It should also be noted that Lense and Thirring in 
Ref. 1 and Soffel in Ref. 4 overestimated the gravitomagnetic precessions of the 
Galilean satellites because they modeled Jupiter as a uniform sphere by assuming 
a = 0.4. 

In order to compute the secular LT shifts ATlt on the Jovian Galilean moons 
it is important to consider that such a system is very complex because of the 
strong mutual perturbations of one satellite on another. In particular, the Laplacian 
resonance 

A^i - 3N2 + 2N3 = (7) 

establishes a relation among the Keplerian mean mean motions^ N of lo, Europa and 
Ganymede. This has consequences also on the LT effect. Indeed, the gravitomagnetic 

''For the physical parameters of Jupiter and its moons see http://ssd.jpl.nasa.gov/sat_gravity.html 
''These ephemerides are available at http://www.imcce.fr. 

°It amounts to 2/3 ~ 0.6 for a hollow spherical shell and to 2/5 = 0.4 for a homogenous sphere. 
''By defintion, the mean mean motion is the secular component in the mean longitude expression. 
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force of Eq. (1) is a small correction with respect to the Newtonian monopole 
and, as usual in perturbation theory, it must be evaluated onto the unperturbed 
Keplerian path for which the velocity entering Eq. (1) is approximately v ~ aN. 
Thus, the gravitomagnctic components of the equations of motion for lo, Europa 
and Ganymede are not independent but they are coupled due to Eq. (7). More, 
secular variations of nodes and pericentres will affect the mean mean motions. To 
preserve the Laplacian relation, variations on mean anomalies are also expected. 
Hence, longitude evolutions of lo, Europa and Ganymede have to be studied. This 
means that differences with respect to the case of independent motions, illustrated 
in Table 2, are expected. This feature has been investigated numerically in the 
following way. 

Table 2. Secular LT shifts, in km, over 125 years 
calculated with Eq. (3) and Eq. (5) by neglect- 
ing the mutual perturbations of one satellite on an- 
other. For the Jupiter's angular momentum the value 
S = 4.33 X 10^* kg m^ s"^ has been adopted. 

lo Europa Ganymede Callisto 

ATlt (kin) -28 -11 -4 -1 



3.2. Numerical simulations 

In order to numerically investigate the effect of the Jovian gravitomagnctic field on 
the orbital motion of the four Galilean satellites, we have performed several numer- 
ical integrations. 

The adopted software is called NOE (Numerical Orbital Elaboration) and is inspired 
from a former work presented in Rcf. 18. It was developed at the Royal Observatory 
of Belgium mainly for natural satellites ephemerides purpose. It is a N-body code 
which incorporates highly sensitive modeling and can generate partial derivatives. 
The latter ones arc needed to fit the initial positions, velocities and other param- 
eters to the observation data. We used the exact modeling, initial conditions and 
parameter values that were used during LI ephemerides elaboration. Let us recall 
that this model introduces the Jovian gravity field by mean of J2, J4, Je coefficients, 
the satellites mutual perturbations (including their respective oblateness coefficients 
J2 and C22) and the Solar perturbation using DE406 ephemerides''. The selection ot 
these perturbations was done after a careful study of the magnitude expected from 
a large set of usually neglected perturbations (planetary and satellites precessions, 
Jovian J3 coefficient ...), but did not considered LT effect at that time. A pertur- 
bation was added in the model only when found significant from the observations 



Hence, it already contains the node and pericenter secular variations. So mean mean motion is 
different from the averaged mean motion (that contains only the mean anomaly). See Ref. 25. 
''In particular, the indirect planetary perturbations are implicitely introduced. 
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accuracy. In particular, secular variations inducing less than hundreds of kilometers 
after one century on the satellite longitudes have generally been considered negligi- 
ble. Indeed, these latter can be easilly absorbed by tiny changes on initial satellite 
positions and velocities. On the other hand, small perturbations may be retained if 
damping enough known frequencies in the system or adding new ones. 

LT effect has been introduced by means of Eq. (1) and Eq. (2). For an explicit 
formulation of all the equations used, we refer to Ref. 18. Initial conditions and 
parameter values are available in Ref. 19. In particular, LI ephemerides result from 
the fit of a high sensitive model to observations covering a time span from 1891 to 
2003. 

The integrator subroutine is the one of Everhart^^, called RA15. It was chosen 
for its speed and accuracy. During the integration a constant step size of At = 0.08 
day was used. To increase the numerical accuracy during the fit procedure (see 
subsection 3.3) we performed integration over ±62.5 years instead of +125 years. 
The numerical accuracy of our simulations is at the level of several meters. 

In Fig. (1) we plot the differences on satellite distances between a numerical 
simulation including LT effect and a second simulation (using the same set of ini- 
tial conditions) but neglecting LT effect. The adopted time span is 125 years and 
covers roughly the modern observation period. In a first step we followed our for- 
mer analytical formulation and so neglected the mutual perturbations (by nullifying 
the satellite masses). In particular, the Laplacian resonance was not present in the 
system at this step. Results are shown in the left panel of Fig. (1). 

As expected, the obtained secular shifts agree with the analytical calculations 
of Table 2. 

In a second step, we reintroduced the satellite masses values (and, thus, the 
Laplacian resonance). Results are shown in the right panel of Fig. (1). As it can 
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Fig. 1. Differences in distance for the four Galilean satellites between a simulation neglecting the 
LT effect and a simulation including the LT effect. The case without the Laplacian resonance is 
shown on the left, while it has been included on the right. The horizontal axes are in years and 
the vertical axes are in km. In both cases the initial conditions and parameter values have been 
taken from Ref. 19. 
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be seen, the shifts remain unchanged for Calhsto at the level of 1.4 km. Rather 
small change appear for Ganymede with a slight decrease from 4.4 km to 4.2 km; 
lo and Europa experience the most important changes passing from 28.6 km to 27 
km and from 11.3 km to 16.3 km, respectively. These differences with our numerical 
estimations are induced by the Laplacian resonance as explained in Section 3.1. 

3.3. Measurability of the Lense-Thirring effects with the current 
ephemerides 

Galilean satellite observations are rather different from those for the artificial satel- 
lites, in the sense that only the right ascension and declination for each natural 
satellite are available. Neither estimation of the Jupiter-satellite distances nor of 
satellite velocities are possible. Hence, reconstruction of observed Keplerian ele- 
ments is not possible. Moreover, the observations have a quite smaller coverage by 
satellite revolutions and are not dispatched equally over the years. This is why the 
method used to detect LT effects with LAGEOS satellites cannot be applied here. 
However, and to estimate the possible detection of LT effect among Galilean system 
in a realistic way, we decided to apply a fit procedure of our numerical simulation 
neglecting LT effect on the simulation including LT effect (acting like some real 
observations) . The induced residuals can then be anticipated as part of the real ob- 
served residuals. If these former are found signicant (at least one order of magnitude 
less), then one can deduce that the effect is detectable. 

This method was already applied successfully to estimate tidal effects among 
the Galilean system in Ref. 

We performed the fit procedure following ephemerides elaboration. As for Ll 
ephemerides, we fitted only the initial positions and velocities of each satellite. A 
sample of 628 daytimes with a 72.8 days time step was used. Differences on Cartesian 
positions for all satellites have been fitted, with no weights assigned. 
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Fig. 2. Differences in distance for the four Galilean satellites between a simulation neglecting the 
LT effect and a simulation including the LT effect. The initial conditions have been fitted. The 
horizontal axes are in years and the vertical axes are in km. 
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In Fig. (2) we show the residuals induced by LT effect after fit and which account 
in ephemerides residuals. It can be noted that a large part of the secular LT signature 
has been removed from the differential time series over 125 years to a level of several 
meters (the level of accuracy of our simulations), while an interesting, long period 
pattern appears for Europa with an amplitude of 20 meters. Most of this variation 
appear on the inclination and the mean longitude (graphs not shown here). 

In Table 3 the corrections delivered by our fit procedure and applied to the initial 
positions and velocities are given. Table 4 presents the same corrections converted 
into Keplerian elements. In the case of Callisto, significant changes appear only 
on O and lj with opposite numerical values. That is consistent with our analytical 
approach. For the three other satellites, significant variations appear both on O and 
oj and A4, which is a consequence of the Laplacian resonance. 

Table 3. Corrections to the initial values of the satellite positions (AU) and velocities (AU/d) 
in a J2000 Earth mean equatorial frame centered on Jupiter. 



Satellite 



lo's position 
lo's velocity 
Europa's position 
Europa's velocity 
Ganymede's position 
Ganymede's velocity 
Callisto's position 
Callisto's velocity 



X 



-1.890118053390E-012 
-4.719160826005E-012 
4.161540903546E-013 
1.710567885254E-012 
1.634014971940E-012 
2.916018121412E-012 
-4.627027927472E-014 
1.145268775648E-013 



y 



-4.300639705467E-013 
-3.184161051147E-012 
7.473039114608E-013 
-1.396029125633E-012 
-3.859797030603E-012 
-1.678258226833E-013 
-5.146803122579E-013 
-4.073677159488E-013 



z 



-1.844960816066E-013 
-1.966170498982E-012 
-1.666243857299E-012 
8.221409664166E-013 
-8.237275878758E-013 
1.368721976036E-012 
4.946080437578E-013 
8.541910668374E-013 



Table 4. Corrections to the initial values of the satellites' Keplerian orbital elements. 



Element 



lo 



Europa 



Ganymede 



Callisto 



a (AU) -7.846973938685E-013 

e 5.625520149536E-010 

i (rad) 2.998455097495E-011 

n (rad) -3.400988557444E-008 

u) (rad) 7.4212940681662E-008 

A (rad) -2.520650355108E-010 



-6.236755903388E-013 

2.509390056859E-010 
-4. 1850009742384E-010 

1.728156151514E-008 
-5. 1427109237778E-009 

5.146869597183E-010 



-4.882674126127E-013 
-1.085141494890E-010 
2.2864355504401E-010 
-1.537507277049E-008 
1.2778481739417E-007 
3.610107768281E-010 



-3.689635402759E-013 
-2.457372673403E-011 
1.5703819668256E-010 

3.433878248415E-008 
-3.471 1791663256E-008 

2.891642481017E-011 



4. Conclusions 

In view of the latest developments in our knowledge of the dynamics of the Jo- 
vian system of Galilean satellites, in this paper we investigated the influence of the 
gravitomagnetic field of Jupiter on their motion and how difficult a detection of 
the Lense-Thirring effect would be still now. We pointed out the presence of sec- 
ular drifts on mean longitudes with a highest amplitude of few tens of kilometers. 
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However, most of this effect will vanish during the fit procedures that are used in 
ephemerides elaboration. Resulting residuals have a maximum magnitude of 20 me- 
ters (over 125 years). This appears negligible with today's best observations which 
have an accuracy of a few tens of kilometers. A spacecraft orbiting Jupiter or the 
adjunction of new more accurate observations of the Galilean satellites, is definitely 
required in order to reveal the Lense-Thirring effect in the Jovian system. 
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